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ABSTRACT 

We present new far-infrared (70 — 500 /im) Herschel PACS and SPIRE imaging ob- 
servations as well as new mid-IR Gemini/T-ReCS imaging (8.7 and 18.3 /im) and 
spectroscopy of the inner Lindblad resonance (ILR) region (R < 2.5 kpc) of the spiral 
galaxy NGC 1365. We complemented these observations with archival Spitzer imaging 
and spectral mapping observations. The ILR region of NGC 1365 contains a Seyfert 1.5 
nucleus and a ring of star formation with an approximate diameter of 2 kpc. The strong 
star formation activity in the ring is resolved by the Herschel/PACS imaging data, as 
well as by the Spitzer 24 /mi continuum emission, [Neil] 12.81 /im line emission, and 6.2 
and 11.3 /im PAH emission. The AGN is the brightest source in the central regions up 
to A ~ 24 /zm, but it becomes increasingly fainter in the far-infrared when compared 
to the emission originating in the infrared clusters (or groups of them) located in 
the ring. We modelled the AGN unresolved infrared emission with the CLUMPY torus 
models and estimated that the AGN contributes only to a small fraction (~ 5%) of the 
infrared emission produced in the inner ~ 5 kpc. We fitted the non-AGN 24 — 500 /im 
spectral energy distribution of the ILR region and found that the dust temperatures 
and mass are similar to those of other nuclear and circumnuclear starburst regions. 
Finally we showed that within the ILR region of NGC 1365 most of the on-going star 
formation activity is taking place in dusty regions as probed by the 24 /im emission. 

Key words: Galaxies: evolution — Galaxies: nuclei — Galaxies: Seyfert — Galaxies: 
structure — Infrared: galaxies — Galaxies: individual: NGC 1365. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
f E-mail: aalonso@ifca.unican.es 



1 INTRODUCTION 

The fueling of the nuclear and circumnuclear activity of 
galaxies has been a topic of extensive discussion. Such activ- 
ity not only includes the accretion of matter onto a super- 
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Table 1. Summary of Herschel photometry observations 



Obsid Instrument Bands Start time Duration 

(urn) (UTC) (s) 



1342201436 


SPIRE-P 


250, 350, 500 


2010-07-14 20:20:45.0 


999 


1342222495 


PACS-P 


70, 160 


2011-06-11 12:59:25.0 


1217 


1342222496 


PACS-P 


70, 160 


2011-06-11 13:20:45.0 


1217 


1342222497 


PACS-P 


100, 160 


2011-06-11 13:42:05.0 


1217 


1342222498 


PACS-P 


100, 160 


2011-06-11 14:03:25.0 


1217 



massive black hole with the accompanying active galactic 
nucleus (AGN), but also the presence of intense nuclear and 
circumnuclear starbursts. Both types of activity require gas 
to be transported from the host galaxy on physical scales 
of a few kiloparsecs down to less than one kiloparsec for 
the nuclear starburst activity and even further in for the 
nuclear activity. Interactions, mergers, and large-scale bars, 
among others, have been proposed as possible mechanisms 
to transport gas from kiloparsec scales t o the nuclear and 
circumnuclear regions (see the review bv l.Togeel liool and 
references therein). 

In isolated galaxies with a large-scale bar the gas is 
believed to flow inwards between corotation and the inner 
Lindblad resonance (ILR). Indirect evidence of this is the 
presence of star formation rings and hot spots near the ILR 
of barred galaxies (see the review of lButa fc Combesi ri996). 
The direction of the flows inside the ILR is generally out- 
wards. This implies the existence of a gravity torque barrier 
at this resonance. However, if there is a "spatially-extended" 
ILR region, this translates into the existence of an inner 
ILR (IILR) and outer ILR (OILR). Numerical simulations 
predict that in this case a gas ring can be formed in be- 
tween these two resonances. Furthermore, in the scenario of 
a double ILR the dynamical decouplin g of an embedded nu- 
clear bar can dri ve the gas further in dShlosman et al.ll 19891 : 
iHunt et al.ll2008l ; ICarcfa-Burillo et aUl2009h . The combined 
action of gravity torques due to embedded structures (bars- 
within-bars) and viscous torques could be a viable mech- 
anism to drive the gas to the inner few parsecs and feed 
the A GN (jGarcfa-Burillo et alJ 120051 : iHopkins fe Quataertl 
l2010h . 

NGC 1365 is a giant isolated barred galaxy at a dis- 
tance of 18.6 Mpc ijLindbladl Il999l . th erefore 1" = 90 pc) . 
This galaxy hosts a Sey fert 1.5 nucleus (Sc hulz et al. II 199*91 ) . 
I Jubgwiert et al. I (|l997t ) showed that in NGC 1365 there is 
also a nuclear ba r (R < 10") embedded in the the large-scale 
bar (R ~ 100'M. lLindbiad" et al] l|l996t ) used hydrodynami- 
cal simulations to reproduce the kinematics and the offset 
dust lanes along the large-scale bar of this galaxy with an 
outer ILR at a radius of J?oilr = 27" = 2.4 kpc. There is 
also an inner ILR at a radius of Rulr ~ 0.3 kpc. Henceforth 
we refer to the ILR region as the region interior to the OILR 
of NGC 1365. As predicted by simulations, there is a ring of 
star formation inside the ILR region of this galaxy. 

The star formation activity in the central regions of 
NGC 1365 has been rev ealed by the presence of hot spots 
JSersic fc Pastorizalll965 | ). intense Ha emission (AD oin et al] 
Il98ll : iKristen et all Il997f ). non-thermal radio continuum 
source s associated with Hn regions and supernov a rem- 
nants jSandavist et al.lll995l : iForbes fc Norris! Il998l ^. large 



amounts of molecular gas (|Sakamoto et al ] |2007l ), and point- 
like and di ffuse extended X-ra y emission not associated with 
the AGN (|Wang et al] 12009 V Moreover, there is evidence 
that a significant fraction of the circumnuclear star forma- 
tion activity might be obscured based on the prominent dust 
lane crossing the nuclear region (see Fig. [1] right panel) and 
the bright and extended mid-infrared (mid-IR) emi ssion in 
this region jTelesco et al.lll993l : lG~alliano et al]|2005h . 

In this paper we present new far-infrared (far-IR) imag- 
ing observations of NGC 1 365 performed with the ESA 
Herschel Space Observatory (|Pilbratt et al] 2010) and new 
mid-IR imaging and spectroscopy obtained with the cam- 
era/spectr ograph Thermal-R egion Camera Spectrograph 
(T-ReCS; iTelesco et al] 1 19981 ) instrument on the Gemini- 
South telescope. The Herschel images were obtained using 
the Photodetector Ar ray Camera and Spectrometer (PACS; 
iPoglitsch et al] |2010T ) and the Spectral and Photometric 
Imaging REceiver (SPIRE; iGriffin et aifalOt ) instruments. 
We also use archival S pitzer data taken with the Infrared 
Array Camera (IRAC; iFazio et al.ll2004l). the Multi-Ban d 
Imaging Photometer for Spitzer (MI PS: Iffieke et alj|2qo4 ). 
and the InfraRed Spectrograph (IRS; iHouck et alj|2004t ) in- 
struments. Using IR observations to study the nuclear and 
circumnuclear activity in the ILR region of NGC 1365 is 
crucial because the central region is crossed by a prominent 
dust lane that obscures from our view in the optical a sig- 
nificant fraction of emission produced there. This paper is 
organized as follows. We describe the observations in Sec- 
tion 2. We analyze the AGN IR emission and the spatially 
resolved IR emission in Sections 3 and 4, respectively, and 
we summarize our conclusions in Section 5. 



2 OBSERVATIONS 

2.1 Herschel/PACS and SPIRE imaging 

We obtained Herschel far-IR imaging observations of 
NGC 1365 using PACS at 70, 100, and 160 /im and SPIRE 
at 250, 350, and 500 /jm. The data are part of the guaran- 
teed time program entitled "Herschel imaging photometry 
of nearby Seyfert galaxies: testing the coexistence of AGN 
and starburst activity and the nature of the dusty torus" 
(PI: M. Sanchez-Portal). 

The PACS observations were carried out using the stan- 
dard scan map mode that takes two concatenated scan line 
maps at 45° and 135° (in array coordinates), at a speed of 
20 arcsec/sec, each one with 26 lines of 9' length and cross- 
scan step of 20". This mode produces a rather homogeneous 
exposure map within a square region of about 7' x 7'. The 
set of maps were duplicated to observe through both the 
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Figure 1. Herschel/PACS 100 fim image (left panel) of NGC 1365. The image is shown in a square root intensity scale. The displayed 
FoV is 36 0" X 360" and m a tches that of the optical image shown in the right panel. The latter reproduces the BVR image shown in 
figure 1 of Elmcgrcc n et all ll2009h (reproduced by permission of the AAS), for an easy comparison of the large-scale optical and far-IR 
morphologies. The square shows a FoV= 60" X 60" as in Fig. [2] The small ellipse on the optical image represents the approximate size 
of the ILR region studied in this work, whereas the large ellipse shows the corotation radius. The horizontal bar represents 10 kpc. For 
both images, orientation is north up, east to the left. 



70 (J.m ("blue") and 100 [im ("green") filters. Therefore the 
galaxy was observed twice through the 160 u.m ( "red" ) filter. 
With the SPIRE photometer we observed the three available 
bands simultaneously using the "large map" mode, with two 
nearly orthogonal scan maps (2 scan lines each), at a scan 
speed of 30arcsec/sec, and three repetitions. The homoge- 
neous exposure area for scientific use is approximately 8' x 8'. 
Table [T] gives the summary of the observations. 

We reduced the data with the Herschel Interactive 
Processing Env ironment (HIPE) v8.0.1 and Scanamorphos 
|Roussel I |20T3 ) vl5. For the PACS instrument, we used 
HIPE and the Calibration Database V32 to build Level 1 
products. These included detecting and flagging bad pix- 
els, converting the ADU readings to flux units (Jy/pixel), 
and adding the pointing information. We did not attempt 
to perform deglitching at this stage to prevent the bright 
AGN nucleus to be affected by the MMT deglitching process. 
The final maps were built from the Level 1 products using 
Scanamorphos, which performs a baseline subtraction, cor- 
rection of the striping effect due to the scan process, removal 
of global and individual pixel drifts, and finally the map 
assembly using all the nominal and cross-direction scans. 
For SPIRE we used the standard (small) map script and 
Calibration Database v8.1. The processing included detec- 
tion of thermistor jumps in the time line, frame deglitch- 
ing, low pass filter correction, conversion of readings to flux 
units (Jy/beam), temperature drift and bolometer time re- 
sponse corrections, and addition of pointing information. 
We built the final maps using the "naive" scan mapper 



task. Colour corrections (for PACS, see lPoglitsch et al 1l2010l ; 
please refer to the Observer's Manual for the SPIRE ones) 
are small for blackbodies at the expected te mperatures (e.g., 
iPerez Garcia fc Rodriguez Espinosal 1200 ll ) and have been 
neglected. More details on the processing of Herschel data 
are given in Sanchez-Portal et al. (in preparation). 

Figure Q] shows the fully reduced PACS 100 /im im- 
age of NGC 1365 together with the optical BVR image 
from lElmegreen et all (120091 ). Fig. [2] shows the PACS im- 
ages together with the SPIRE 250 /an with a field of view 
(FoV) covering the approximate extent of the ILR region 
of NGC 1365 (see also Section f2 . 5 f) . We performed aperture 
photometry on all the Herschel images using two different 
radii, r = 15" and r = 30". The latter encompasses the ILR 
region, whereas the former includes mostly the ring of star 
formation. Table lists the measured flux densities. 



2.2 Gemini/T-ReCS imaging and spectroscopy 

We obtained mid-IR imaging of NGC 1365 using T-ReCS 
on the Gemini-South telescope on September 8 and 9, 2011 
as part of proposal GS-2011B-Q-20 (PI: N. Levenson). We 
used the Si-2 (A c = 8.74 (im) and the Qa (A c = 18.3 /jm) fil- 
ters and mid-IR standard observation techniques. The plate 
scale of the T-ReCS imaging observations is 0.089"/pixel 
with a FoV of 28.5" x 21.4". The total integration times 
(on-source) were 145 s and 521s in the Si-2 and Qa filters, 
respectively. The Qa filter observations were split between 
the two nights, whereas those in the Si-2 filter were done on 
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arcsec 


-20 






Figure 2. IR view of the approximate extent of the ILR region of NGC 1365. We show the Spitzer/IRAC images at 3.6, 4.5, 5.8, and 
8/^m, the Spitzer/MIPS image at 24 pm, the Herschel/PACS 70, 100, and 160 /im images and the Herschel/SPIRE image at 250 £tm. 
We mark as filled stars the p ositions of the brighte st mid-IR clusters (M2, M3, M4, M5, M6, see text for more details and also Fig. [4} 
dete cted from the grou nd by iGalliano et all i2005T l, the position of the AGN as a filled dot, and the position of the L4 Ha hot spot 
from lAlloin et al.l il98lT) as a filled square (see Table l4l . Orientation is north up, east to the left. The images are shown in a square root 
intensity scale. 



the second night. We observed standard stars immediately 
before or after the science observations in the same filters, to 
both flux-calibrate the galaxy observations and to estimate 
the unresolved nuclear emission (see below). The observa- 
tions were diffraction limited, with a full width half maxi- 
mum (FWHM) of 0.34" in the Si-2 filter and 0.55 - 0.58" 
in the Qa filter, as measured from the standard star obser- 
vations. We reduced the imaging data using the CanariCam 
data reduction package RedCan (Gonzalez-Martin et al., in 



preparation) . We refer the inte rested reader to this work and 
I Ramos Almeida et alj (|2011bl ) for details. Fig. [3] shows the 
fully reduced T-ReCS Qa image resulting from the combi- 
nation of the data taken during the two observing nights. 
The T-ReCS 8.7 (im image (not shown here) shows a similar 
morphology. 

We also retrieved archival T-ReCS spectroscopic obser- 
vations in the iV-band (~ 8 — 13 /im) using a 0.35" slit width 
as part of proposal GS-2009B-Q-19 (PI: M. Pastoriza). The 
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Table 2. Mid and far-IR aperture photometry of the central region of 
NGC 1365 



Instrument 


Ac 


Pixel size 


FWHM 


/, 0=15") 


f v (r=30") 




w 


(arcsec) 


(arcsec) 


(Jy) 


(Jy) 


MIPS 


24 


2.45 


5.9 


7.2 


8.7 


PACS 


70 


1.4 


5.6 


85.5 


102.7 


PACS 


100 


1.7 


6.8 


110.1 


141.7 


PACS 


160 


2.8 


11.3 


87.4 


123.7 


SPIRE 


250 


6.0 


18.1 


33.4 


53.3 


SPIRE 


350 


10.0 


25.2 


10.4 


21.0 


SPIRE 


500 


14.0 


36.9 




5.6 



Note. — The reported values of the FWHM are the nominal values. The errors 
in the aperture photometry are dominated by the photometric calibration 
uncertainty of the instruments that is typically 10%. 



total on-source integration time was 600 s. We reduced the 
galaxy and corresponding standard star observations us- 
ing the RedCan package fo llowing the steps described in 
lAlonso-Herrero et alj l|201ll ). Finally we extracted the nu- 
clear spectrum as a point source. 

To estimate the nuclear unresolved emission from 
the mid-IR imaging data we followed the point spread 
function (PSF) sca li ng t ec hnique implemented by 
iRamos Almeida et alj l|2009l . l2011bl ). This unresolved 
emission which is assumed to represent the torus emission 
(see Section 13. To do so, we scaled the observation of 
the corresponding standard star to the peak of the nuclear 
emission of the galaxy in each of the two filters. This 
represents the maximum contribution of the unresolved 
source (100%), whereas the residual of the total emission 
minus the scaled PSF accounts for the extended emission. 
In both filters we found that a 90% PSF scaling provided a 
realistic estimate of the extended emission. The estimated 
errors in the T-ReCS unresolved flux densities reported 
in Table [3] are 15% and 25% in the Si-2 and Qa filters, 
respectively, and account for both t he flux calibration and 
PSF subtraction uncertainties (see IRamos Almeida et al.l 
120091 . for more details). The unresolved 8.7 fim emission 
computed this way is in good agreement with the flux 
density at the same wavelength obtained from the T-ReCS 
nuclear spectrum. 



2.3 Archival Spitzer/IRAC and MIPS imaging 

We retrieved from the Spitzer archive imaging data of 
NGC 1365 observed with all four IRAC bands (3.6, 4.5, 
5.8, and 8/xm) and with MIPS at 24 ^m (Program ID: 
3672, PI: J. Mazzarella) . These observations were part of 
Th e Great Observator ies All-Sky LIRG Survey (GOALS, 
sec Armus et al.l [2009). We retrieved the basic calibrated 
data (BCD) from the Spitzer archive. The BCD process- 
ing includes corrections for the instrumental response (pixel 
response linearization, etc.), flagging of cosmic rays and sat- 
urated pixels, dark and flat fielding corrections, and flux 
calibration based on standard stars. We combined the BCD 
images into mosaics using the MOsaicker and Point source 
Extractor (MOPEX) software provided by the Spitzer Sci- 
ence Center using the standard parameters. The final mo- 
saics were repixeled to half of the original pixel size of the 
images, that is, the IRAC mosaics have 0.6"/pixel, whereas 



the MIPS 24 /im mosaic has 1.225" /pixel. In Fig.fjwe show 
the Spitzer images with a FoV covering the approximate ex- 
tent of the ILR region, as done with the PACS images and 
the SPIRE 250 /im image. The angular resolutions of the 
IRAC images are between 1.7 and 2" (FWHM) and that of 
the MIPS 24 /im is 5.9". 



2.4 Archival optical ground-based imaging 

We retrieved from the ESO archive optical images obtained 
with the WFI instrument on the MPG/ESO 2.2m tele- 
scope using the narrow- band Halpha/7 filter (A c = 6588.3A, 
FWHM=70A) obtained as part of proposal 065.N-0076 (PI: 
F. Bresolin). We combined a total of 6 images, each of 
them with a 350 s exposure. The plate scale of the images 
is 0.238"/pixel. The filter contains the Ha and [Nil] emis- 
sion lines plus adjacent continuum. Since we use this image 
for morphological purposes only, we did not attempt to ei- 
ther subtract the continuum or calibrate it pho tometrically. 
The p ositions of the Ha hot spots identified bv lAlloin et al.l 
l|l98ll ) in the central region of NGC 1365 are displayed in 
Fig. [3] (upper panel). 



2.5 Alignment of the images 

The alignment of the Spitzer/IRAC and the Gemini/T- 
ReCS images is straightfor ward because the AG N and the 
mid-IR clusters detected bv lGalliano et afl ((2005) are clearly 
identified in all these images. The AGN is also bright in 
the optical image and therefore we used it as our reference. 
In Fig. f3] we present a close-up of the ILR region in the 
IRAC 8 fj,m band and the optical narrow-band Ha image. 
We marked the positions of the bright (designated as M4, 
M5, and M6) and faint mid-IR (designated as M2, M3, M7, 
and M8) clust ers and the nucleus using the relative posi- 
tions given bv lGalliano et all ((2005). We also indicated the 
positions of radio sources and Ha sources (see Section [4. ip . 
Although the Spitzer/MIPS 24 /im image has a poorer an- 
gular resolution when compared to that of IRAC, the AGN 
is still sufficiently bright at 24 /im that can be distinguished 
from the mid-IR clusters. 

The alignment of the Herschel/PACS images is not as 
simple because the AGN does not appear to be a bright 
source in the far-IR. We used the astrometry information 
in the Herschel image headers and the optical position of 
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10 -10 

arcsec 

Figure 3. Gemini/T-ReCS imago of the central region of 
NGC 1365 in the Qa (A c = 18.3 /im) filter. The image has been 
rotated to the usual orientation of north up, east to the left, and 
smoothed with a Gaussian function. We mark the positions of 
the AGN (fil led dot) as well as the mid-IR clusters (star symbols) 
identified bv lGalliano et al] d2005lV 

the nucleus of NGC 1365 given by ISandqvist et all l|l995h : 
RA(J2000)=03 h 33 m 36.37 s and Dec(J2000)=-36°08'25.5", 
for the initial alignment. These coordinates placed the AGN 
to the southwest of the bright source detected in the PACS 
images and in the MIPS 24 ^m image that appears to coin- 
cide with the region containing clusters M4, M5, and M6. 
Additionally we compared the morphologies of the PACS 
70 /im and the MIPS 24 ^m images as they have similar an- 
gular resolutions (see Table [5} ■ We used the mid-IR source 
located ~ 17" N, 16" E of the AGN identified in the IRAC 
and MIPS images and also seen in the PACS 70 and 
100 (im images for a finer alignment. This source appears to 
be coincident with the L4 H II region or Ha hot spot (see 
Fig. Hand Tableg} identified bv lAlloin etafl (|l98ll ) and also 
seen in our archival Ha image. The PACS 160 /an image was 
aligned relative to the other PACS images with the astrome- 
try information in the headers. Finally since the bright mid- 
IR clusters cannot be resolved in the SPIRE 250 /im image, 
we placed the center of image at the position of the AGN. 
We note that the alignment of the PACS images in Figs. [5] 
and |3] is only good to within 1 pixel in each band. Since the 
main goal of this work is to study the processes giving rise 
to the IR emission within the ILR region of NGC 1365, in 
Fig. [5] we show the aligned IR images with a FoV covering 
the approximate extent of this region. We do not show the 
SPIRE 350 and 500 /im images because of the small number 
of pixels covering the ILR region of NGC 1365. 

2.6 Spitzer/IRS spectral mapping 

We retrieved from the Spitzer archive low spectral resolu- 
tion (R ~ 60 - 126) observations (Program ID: 3269, PI: J. 
Gallimore) of NGC 1365 obtained with the spectral map- 
ping capability of IRS. These observations were part of the 
Spitzer study of the spectral energy distributions (SEP) o f 
the 12 /im sample of active galaxies (iGallimore et al.ll2010l ). 




arcsec 

o n — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — r 

C\2 




arcsec 

o n — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — r 

C\2 




arcsec 

Figure 4. Close-ups of the inner 36" X 36" region showing the 
location of the AGN ( filled dot), the mid-I R clusters M2...M8 
(filled star symb ols) of Galliano et al. (20Q5|), the radio sources 
(open circles) of ISand qvist et"a l I lll995l ), and the Ha hot spot s 
(filled squares) of lAlloin et al] jl98ll Tand iKristen et al] jl997l) . 
The upper panel is the optical ESO Halpha/7 narrow-band image, 
the middle panel is the IRAC 8 fim image, and the lower panel 
the PACS 70 fim image. The beam of the images is shown on the 
lower left corner of each image, approximated as the FWHM of a 
Gaussian function. ^ 000Q RAg ^ MNRAS 000 000-000 
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Table 3. AGN emission 



Wavelength /„ Method 

(/im) (mjy) 



8.7 


203 ± 30 


Imaging (unresolved) 


13.0 


400 ± 60 


Spectroscopy 


18.3 


818 ± 205 


Imaging (unresolved) 


24 


1ZOO -500 


BC fit 


70 


7 o 4 +1482 
' °*-422 


BC fit 


100 


971 +632 
z,i -163 


BC fit 


160 


< 78 


BC fit 



References. — The quoted uncertainties for the BC fit 
fluxes are the ±lo" uncertainties, as discussed in Sec- 
tions [XT] and GLU 



PAH [Aril] PAH PAH [ArIlI]H 2 [SIV] PAH PAH [Nell] 













I 1 - 














-d 










j \J AGN (x2.8) ; 


0J 












N 
























s 












u 
o 












S5 








msSy? M2+M3 l 










\ M4 












-If M5+M6 











5 10 



Rest wavelength (/im) 

Figure 5. Spitzer/IRS SL1+SL2 5 — 15 /an spectra of selected 
regions (see Table [SJ normalized at 13 /im. The spectrum of the 
AGN has been shifted up for clarity. The most prominent spectral 
features are marked. 

The observations were obtained with the Short-Low (SL1; 
7.5 — 14.3 /mi and SL2; 5.1 — 7.6 /mi) and Long-Low (LL1; 
19.9 - 39.9 /im and LL2; 13.9 - 21.3 /mi) modules. The 
plate scales of the SL and LL modules are 1.8"/pixel and 
5. l"/pixel, respectively. 

The data were processed using the Spitzer IRS pipeline. 
The IRS data cubes were assembled using CUBISM (the 



Table 4. Summary of sources in the circumnuclear region 



Name 


Spectral 


Rel. Position 


Ref. 




Range 


arcsec, arcsec 




M2, L2 


mid-IR, Ha 


-4.7, -5.1 


1, 2 


M3, L3 


mid-IR, Ha 


-5.4, -2.6 


1, 2, 3 


M4, D 


mid-IR, radio 


0.4, 7.1 


1, 4 


M5, E, L12 


mid-IR, radio, Ha 


2.8, 10.0 


1, 4, 3 


M6, G 


mid-IR, radio 


4.8, 7.0 


1, 4 


M7, H 


mid-IR, radio 


* 


1, 4 


M8, LI 


mid-IR, Ha 


* 


1, 2 


Lll 


Ha 


10, 15 


3 


L4 


Ho 


17, 16 


3 


A 


radio 


-4.1, -2.4 


4 



The positions are relative to that of the AGN and correspond to 
those of the first listed source. *The positions shown in Figs. 3 
and 4 are estimated from the mid-IR images of Galliano et al. 
(2005). References.— 1. Galliano et al. (2005). 2. Kristen et al. 
(1997). 3. Alloin et al. (1981). 4. Sandqvist et al. (1995). 



CUbe Builder for IRS Spectra Maps, ISmith et al l l2007al ) 
vl.7 from the individual BCD spectral images obtained from 
the Spitzer archive. CUBISM also provides error data cubes 
built by standard error propagation, using, for the input un- 
certainty, the BCD-level uncertainty estimates produced by 
the IRS pipeline from deviations of the fitted ramp slope fits 
for each pixel. We used these uncertainties to provide error 
estimates for t he extracted spectra , and the line and contin- 
uum maps (see lSmith et al.] [2007a, for full details) discussed 
in the next two sections. 



2.6.1 Extraction of the ID spectra 

We used CUBISM to extract spectra of regions of inter- 
est using small apertures taking advantage of the angu- 
lar resolution of the spectral mapping observations ob- 
tained with the SL1+SL2 m odules (~ 4" FWHM, see 
iPereira-Santaella et al.ll20T0bl ). We used square or rectangu- 
lar apertures in the original orientation of the SL data cubes 
with sizes of two or three pixels (see Table [5] for the extrac- 
tion apertures used for each region). The selected regions 
include the AGN, and the regions containing the M2+M3, 
M4, and M5+M6 mid-IR clusters (see F ig. 0] for the posi- 
tions) identified bv lGalliano et all (|2005l L We note that we 
did not attempt to apply a point source correction to the SL 
3.7" x 3.7" spectrum of the AGN because we are mostly in- 
terested in the extended features, that is, the polycyclic aro- 
matic hydrocarbon (PAH) features and the [Ne n] 12.81 /mi 
fine structure line. Fig. [S]shows the SL1+SL2 spectra of the 
selected regions normalized at 13 /im. Finally we extracted 
the integrated spectrum of the region covered by the obser- 
vations, i.e. the central 27.8" x 24.0". 

We measured the fluxes and the equivalent widths (EW) 
of the [Nen]12.81 /mi emission line and the 6.2 and 11.3 /im 
PAH features fitting Gaussian profiles to the lines and lines 
to the local continuum. Our measurements of the [Ne n] flux 
for clusters M4, M5, M6 are in go od agreement with those 
reported bv lGalliano et al] (|2008al ) from ground-based high 
angular resolution observations. Since the PAH features are 
broad, it has been noted in the literature that Gaussian pro- 
files might not be appropriate to measure their flux because 
a large fraction of the energy in these bands in radiated in 
the wings. A Lorentzian profile mi ght be a better approx- 
imation (see I Galliano et al.l l2008bj ) to measure their flux, 
and therefore we repeated the line fits with this method. 
Table [5] lists the measurements for the extracted spectra for 
the Gaussian profiles. To illustrate the effect of using differ- 
ent profiles, in this table we give the measured 6.2 to 11.3 /im 
PAH ratio for the two profiles and the selected regions. 



2.6.2 Spectral Maps 

We used CUBISM to construct spectral maps of the most 
prominent features in the SL data cubes, namely, the 6.2 /im 
and 11.3 /im PAH features, and the [Nen]12.81 /im fine 
structu re line. The technique used h ere was very similar to 
that of I Alonso-Herrero et alj l|2009t ) and involves integrat- 
ing the line flux over the user-defined emission line regions. 
Note that, unlike the line measurements in the previous sec- 
tion, the features are not actually fitted with Gaussian or 
Lorentzian profiles, and therefore these maps are only used 
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Table 5. Measurements from the Spitzer/IRS SL1+SL2 spectra 



Region 


Extraction 
Aperture 


6.2 pm PAH feature 
flux EW 


11.3 /im PAH feature 
flux EW 


[Nc n] 12.81 fan 
flux 




PAH6.2 
PAH11.3 

G L 


AGN 


3.7" X 3.7" 


8.6 


0.09 


5.3 


0.13 


1.7 


-0.10 


1.6 1.2 


M2+M3 


5.6" x 3.7" 


21.0 


0.50 


12.7 


0.61 


1.6 


-0.24 


1.7 1.1 


M4 


3.7" x 3.7" 


10.7 


0.45 


5.8 


0.63 


2.3 


-0.85 


1.8 1.7 


M5+M6 


5.6" x 3.7" 


20.0 


0.45 


10.2 


0.40 


6.5 


-0.81 


2.0 1.5 


Integrated 


27.8" x 24.0" 


265.0 


0.42 


161.0 


0.49 


60.5 


-0.45 


1.6 1.3 



Notes. — Fluxes (observed, not corrected for extinction) are in units of xl0~ 13 erg cm" 2 s~ 1 and EW in units of fim 
for measurements done with Gaussian profiles. The typical errors are 10% for the fluxes and 0.05 fim for the EW. 
Sgi is the strength of the silicate feature (see Section 12.6.21 for the definition). The ratio of the 6.2 fim PAH flux to 
the 11.3 /an PAH feature flux is given for fits to the features done with Gaussians (G) and Lorentzian (L) profiles 
(see Section l2.6.1l for more details). 
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10 -10 10 -10 



arcsec arcsec 

Figure 6. Spitzer/IRS observed (not corrected for extinction) 
spectral maps built with CUBISM from the SL data cubes. The 
maps of the 6.2 fim PAH feature (lower left), 11.3 /im PAH feature 
(upper right), and [Ne II] 12.81 fim (lower right) are in units of 
10 — 7 Wm — 2 sr _1 , and the 5.5 fim continuum map (upper left) 
in units of MJy sr — 1 . Only pixels with relative errors of less than 
20% are displayed. The position of the AGN (filled dot) was made 
to coincide with the peak of the 5.5 fim continuum emission. The 
positions of th e M2.. .M6 mid-IR clusters (filled star symbols) of 
iGalliano et"afl J2005T) are also plotted. Orientation is north up, 
east to the left. The apparently faint 11.3 PAH emission in the 
region of clusters M4, M5, and M6 is due to foreground extinction 
(see also the left panel of Fig. 7, and Sections 14.11 and 14. 21 . 



arcsec arcsec 

Figure 7. Spitzer/IRS spectral maps of the strength of the sili- 
cate feature Ssi (left panel) and the observed (not corrected for 
extinction) ratio of the 6.2 fim to the 11.3 fim PAH feature emis- 
sion (right panel). Negative values of Ssi indicate that the feature 
is observed in absorption. Orientation is north up, east to the left. 
Symbols are as in Fig. [6] 

IPereira-Santaella et al] l)2010bh . Briefly, it involves fitting 
the silicate feature from ID spectra extracted in 2pixel x 
2pixel boxes from the IRS SL data cubes. The map is then 
constructed by moving by 1 pixel in the x and y directions 
until the FoV of the IRS SL data cubes is completely cov- 
ered. We measured the apparent str ength of the silicat e 
feature in the IRS spectra following ISpoon et al] (|2007t ): 
Ssi = hi /„(obs)//„(cont), where /„(obs) is the flux density 
at the feature and /„ (cont) is the flux density of the underly- 
ing continuum. The latter was fitted as a power law between 
5.5 and 14 pm. We evaluated the strength of the silicate fea- 
ture at 10 /im. When the silicate strength is negative, the 
silicate feature is in absorption, whereas a positive silicate 
strength indicates that the feature is seen in emission. The 
uncertainties of the measured strengths in the spectral map 
and the extracted ID spectra of previous section are ±0.05. 



for morphological purposes. Since CUBISM does not fit or 
deblend emission lines, the SL [Nell] 12.81 fim map includes 
some contribution from the 12.7 fim PAH feature. We also 
built a continuum map at 5.5 fim. The Spitzer/IRS SL spec- 
tral maps shown in Fig. [6] were trimmed and rotated to the 
usual orientation of north up, east to the left. We used the 
associated uncertainty maps produced by CUBISM to com- 
pute the relative errors of the spectral maps. 

Finally we constructed the map of the silicate fea- 
ture, which is shown in Fig. [7] following the technique of 



3 AGN INFRARED EMISSION 

A number of works have studied in detail the mid-IR emis- 
sion of the AGN hosted by NGC 1365 using high angu- 
lar resolution observations, including imaging, spectroscopy, 
and interferometry. Our T-ReCS nuclear (~ 0.35") mid- 
IR spectrum of NGC 1365 (see Fig. [SJ is an almost fea- 
tureless continuum with no clear evidence of the presence 
of the 9.7 fim silicate feature or PAH features (see also 
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Table 6. Parameters of the CLUMPY Torus Models 



Parameter 



Interval 



Torus radial thickness Y [5, 100] 

Torus angular width utorus (dog) [15, 70] 

Number clouds along an equatorial ray No [1, 15] 

Index of the radial density profile q [0, 3] 

Viewing angle i (deg) [0, 90] 

Optical depth per single cloud ry [5, 150] 



ISiebenmorgen et all 12004 : iTristram et alj 120091 ). By con- 
trast, the spectrum integrated over several arcseconds is 
rich in spectral features related to star formation activ- 
ity, such as, PAH fea tures and the [Neil] line (see Fig. [5}. 
ITristram et aD l|2009l ) modelled mid-IR interferometric ob- 
servations of the nuclear region and derived a size (diam- 
eter) of 1.1 — 2.7 pc for the 12 /an emitter. Up until now, 
there has been no estimates of the far-IR emission arising 
from the AGN of NGC 1365, and thus the main goal of this 
section is to provide such an estimate. 

On scales of a few arcseconds the AGN of NGC 1365 
is surrounded by a number of mid-IR bright clusters and 
extended and diffuse emission (see Figs. [2] and [4] and Sec- 
tion Moreover, in the far-IR AGN emission becomes less 
luminous for increasing wavelengths when compared to that 
of the clusters (see Fig. [2}. This together with the limited 
angular resolution of the PACS images (i.e., the best angu- 
lar resolution is FWHM~ 5.6" at 70 /mi) prevents us from 
doing aperture photometry on the PACS images to derive 
the far-IR flux densities of the AGN. Instead, in this section 
we estimate the AGN emission by both fitting the nuclear 
near- and mid-IR SED with torus models and doing a spec- 
tral decomposition of the Spitzer/IRS nuclear spectrum. 



3.1 Fit to the nuclear SED and mid-IR spectrum 
using clumpy torus models 

i Honig 



Clumpy torus mo dels (e.g., IHonig: et al. 20061: 

ISchartrnann et aT1l2008l ) and in particular the lNenkova et al.l 
[|2008al lbh models (also known as clumpy models) have 



been shown to fit well the near a nd mid-IR SEDs and 
spect r oscopy o f Seyfert ga l axies dRamos Almeida et al.l 
20091. l2011bl: iMason et al'l 120091: iNikutta et all 120091: 



Alonso-Herrero et all l201ll : ISales et al.l l201ll : iLira et all 
20121 ). The clumpy torus models are described by six 



parameters to account for the torus geometry and the 
properties of the dusty clouds (see Table [6]|. The geometry 
of the torus is defined with the torus radial thickness Y, 
which is the ratio between the outer and inner radii of 
the torus Y — R /rJ^\, and the width of the angular 
distribution of the clouds atoms, that is, the angular size of 
the torus as measured from its equator. The clouds have an 
optical depth tv and are arranged with a radial distribution 

1 The inner radius of the torus in these models is set by the dust 
sublimation temperature, which is assumed to be T su b ~ 1500 K, 
and the AGN bolometric luminosity L DO i(AGN). Then the dust 
sublimation radius in pc is R d = 0.4L bol (AGN) - 5 , where the 
AGN bolometric luminosity is in units of 10 45 ergs -1 . 



3 0.1 r 




Rest wavelength (^in) 

Figure 8. Best fit to the nuclear SED (blue dots) and mid-IR 
spectroscopy (blue line) of NGC 1365 (MAP and median fits are 
the solid and dotted lines, respectively) using the CLUMPY torus 
models (including an AGN component). The shaded region repre- 
sents the range of acceptable models with a 68% confidence level. 
The inset shows the spectral region around the 9.7 /im silicate 
feature. 



expressed as a declining power law with index q (oc r q ), 
with A^o giving the mean number of clouds along a radial 
equatorial ray. The last parameter is the viewing angle to 
the torus i. We refer the reader to the sketch of the cl umpy 
model geometry in Figure 1 oflNenkova et all (l2008ah. 

We updated the fit of I Ramos Almeida et al. ( 2011bl ) 
with two main differences. First, we used the higher angular 
resolution mid-IR unresolved flux densities inferred in Sec- 
tion [2]2] and also included the T-ReCS mid-IR spectrum for 
the fit. Second, this galaxy is clas sified as a Seyfert 1 .5 (it 
shows a broad component in H/3, ISchulz et al. |[l999l ) and 
therefore there is an unobscured view to the AGN. To ac- 
count for it, we added an AGN component to the resulting 
torus model (see below) and allowed for a s mall amount 
of fo reground extinction (Av < 5 mag, see lAlloin et al.l 
1981), to redden the fitted models. We also included in 
the nuclear SEP th e HST/NICMOS 1.6 pm unresolved flux 
jCarollo et alj|2002h . 

Si milarly to our previous wor k (|Ramos Almeida et al.l 

l2011bl : rAlonso-Herrero et al.ll201ll ). we took a Bayesian ap- 
proach to both fit the torus models to the data and de- 
rive meaningful confidence levels for the fitted parame- 
ters. To this end, we used the BayesClumpy (BC) fit- 
ting routine, an interpolated version of th e clumpy models 
|Asensio Ramos fc Ramos Almeidal 12009). The new version 
of BC allows to fit both photometric points and spectra. 
For the modelling we used the whole range of parameters 
probed by the clumpy models. We list them in Table [6] 
Unlike in some of our previous work, here we allowed for the 
torus extent to vary in the range Y = 5 — 100, instead of re- 
stricting the models to small tori. This is because larger tori 
may be required to reproduc e the far-IR emission of AGN 
(R amos Almeida et alj|2011al ). In addition to the six torus 
model parameters to be fitted plus the foreground extinc- 
tion, there is an extra parameter to account for the vertical 
displacement needed to match the fluxes of a given model 
to the observations. This vertical shift, which we allow to 
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vary freely, scales wit h the AGN bolometric luminosity (see 
iNenkova et~ai1l2008bh . We assumed uniform distributions in 
the ranges given in Table HI] for the prior distributions of the 
torus model parameters. 

The BC fitting routine translates the probability distri- 
butions of the fitted torus model parameters into two model 
spectra. The first corresponds to the maximum-a-posteriori 
(MAP) values that represent the best fit to the data, while 
the second is the model produced with the median value of 
the probability distribution of each parameter. Fig. [8] shows 
the best fit to the nuclear near and mid-IR SED and the 
mid-IR spectrum of the NGC 1365 AGN. We also plot in 
this figure the range of acceptable torus models at the 68% 
confidence level. As can be seen from this figure, the fit to 
the photometric and spectroscopic data is very good. Specif- 
ically, the clumpy models rep roduce well the flat s ilicate 
feature at around 9.7 fim (see INenkova et all l2008bl . their 
figures 16 and 17). 

We can now compare the inferred torus model param- 
eters with observations. We start with the AGN bolometric 
luminosity. As explained above, we chose to keep the scaling 
of the torus models to the observations as a free parame- 
ter. Because the scaling is proportional to the AGN bolo- 
metric luminosity, when compared with the observations 
it provides information ab out the goodness of the fit (see 
lAlonso-Herrero et al.|[201ll , for more details). The derived 
AGN bolometric lumin osity from our fit is Lboi(AGN) = 
2.6 ± 0.5 x 10 43 ergs~\ lmsaliti et all (|2005l ) infered a 2 - 
10 keV luminosity in the range 0.8 - 1.4 x 10 42 ergs _1 (cor- 
rected for the assumed distance in this work). Our derived 
AGN bolometric luminosity is consistent with the values 
from X-rays after applying a bolometric correction. 

We can also compare the torus radius R a = 5 jli' 5 pc and 
torus angular width atoms — 36^g 4 deg from the modelling 
with observations. The width of the angular distribution of 
the clouds in the torus is related to the opening angle of the 
ionization cones, with some de pendence with the geometr y 
of the gas (see discussion in I Alonso-Herrero et al. 1 120111 ). 
The half op ening ang l e of t he NGC 1365 ionization cone 
modeled by ILindbladl l|l999l ) is O CO nc = 50deg (measured 
from the pole), which is compatible with the derived angular 
size of the torus. The size (diameter) of the 12 fj,m emitting 
source inferred from the modelling of the mid-IR interfer- 
ometric observations of this galaxy is si2 M m = 1-1 — 2.7pc 
l|Tristram et al.l 120091 ). We note that the size of the 12 /im 
emitting source, which traces the warm dust within the 
torus, is expected to be smaller than the true size of the 
torus. This is because the 12 fim emission is rather insensi- 
tive to cooler material further from the nucleus (see discus- 
sion in lAlonso-Herrero et al.l l201ll ) . Given the good agree- 
ment with the observations, we conclude that the fitted torus 
model provides a good representation of the AGN torus 



3.2 AGN torus emission in the far-IR 

To estimate the AGN far-IR flux density in the PACS bands, 
we extrapolated the fitted torus model beyond ~ 20 fim. As 
can be seen from Fig. [8] the torus model flux density peaks 
at 40—50 fim and falls off rapidly at longer wavelengths. This 
is sim ilar to results for other AGN (see e.g.. lMullanev et al. I 
l201ll ). In Table |3] we list the extrapolated AGN far-IR flux 



densities together with those in the mid-IR from the T-ReCS 
imaging and spectroscopy. The ranges given in this table 
for the far-IR fluxes from the BC fit of the NGC 1365 take 
into account the ±ler confidence region of the fitted clumpy 
torus models. 

A quick comparison between the extrapolated AGN far- 
IR fluxes and the aperture photometry of the central re- 
gions (Table [2j clearly shows that the AGN contribution at 
these wavelengths is small and decreases for increasing wave- 
lengths. Within the central 5.4 kpc (r = 30") the AGN is 
predicted to contribute 15% (+8%, —6%) of the total 24 )im 
emission. At 70 /im the AGN accounts for at most 2% of the 
emission in the same region and less than 1% at 100 and 
160 (im. Using the flux den sities measured for r = 30" and 
the following equation from lDale fc Heloul (|2002l ). 

£3-iooo M m = 1.559 i/L„(24/«n) + 0.7686 uL u (70fjm) + 

1.347 vL v (160fj,m) (1) 

we estimate a 3 — 1000 /im luminosity within the ILR re- 
gion of ~ 9 x 10 10 L Q . The ratio between the AGN bolo- 
metric luminosity and the IR luminosity of this region is 
then 0.05 ± 0.01. This is in good agreement with find- 
ings for the majority of local LIRGs hosting an AGN 



mgs tor tne majority ot k> 
J Alonso-Herrero et alBoil ). 



4 SPATIALLY-RESOLVED IR EMISSION IN 
THE ILR REGION 

4.1 Morphology 

The central region of NGC 1365 is a bright source of 
mid- and far-IR emission (see e.g., Telesco et alj Il993l ; 



iGalliano etlll 120051 ; IWiebe et ail 120091 ). although its con- 
tribution to the integrated emission of the galaxy decreases 
significantly at longer wavelengths. The mid-IR morphology, 
especially in the IRAC 5.8 and 8 /jm and in the T-ReCS im- 
ages, resembles that of a circumnuclear ring of star forma- 
tion with a diameter of approximately 1.8 kpc (see Figs. [2] 
and [3]). The AGN and the bright mid-IR clusters M4 M5, 
and M6 identified (see TableiJ bv lGalliano et al.l (|2005l ) out 
to A ~ 12 /im are also clearly identified in the Gemini/T- 
ReCS image at 18.3 /im (see Fig. [3}- The region to the 
southwest of the AGN is more diffuse and similar to the Ha 
morphology (see Fig. [4] upper panel). We do not, however, 
identify any bright compact sources there. This diffuse emis- 
sion is reminiscent of the optical morphology of the region 
containing superstar clusters SSC3 and SSC6 along with a 
number of less luminous optical clusters (|Kristen et al.lll997l ; 
ILindbladl Il999l ). 

In the Spitzer/MIPS 24 fim image the emission from the 
AGN is still differentiated from that of the mid-IR clusters, 
although the latter are not resolved individually because of 
the decreased angular resolution. The Herschel observations 
of NGC 1365 show that the far-IR AGN emission becomes 
faint when compared to that of the clusters. This is not only 
because the clusters in the central region of NGC 1365 are 
still intrinsically bright in the far-IR, but also because in gen- 
eral the AG N emission is observed to fall steeply in th e far- 
IR (see e.g., iNetzer et aT]|2007i ; iMullanev et al. Il201ll . and 
references therein) . Moreover, in general the far-IR emission 
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of A GN appears to be dominate d by a starburst compo- 
nent ( Hat ziminaoglou et al.ll2010h . It is also apparent from 
Fig. [5]that the group of clusters M5...M8 are brighter than 
the group of M2 and M3 in the far-IR, as is also the case 
in the mid-IR. This implies, as will be discussed further in 
Section [4.51 that the obscured star formation rate (SFR) is 
higher in the former group than in the latter. Moreover, 
the morphology of the SFR tracers (e.g., the 8, 24, and 
70 /im emissions) is similar to that of the mol ecular gas, espe- 
cially that of the rare C ls O(J=2-l) isotope l|Sakamoto et al.l 
120071 ). where a "twin-peak" morphology is observed. Such 
molecular gas twin peaks are observed in barred galaxies 
|Kennev et all 1 1992ft and are formed where dust lanes in- 
tersect nuclear rings of star formation, as is the case of 
NGC 1365. 

We can also compare the IR morphology, in particular 
at 8, 24, and 70 /im, of the ILR region with that o f Ho[f|. Since 
all th e se emissions probe the SFR in galaxies llKennicutt I 



19981; iKennicutt et all 120091; IXionso-Herrero et all |20o5 



Calzetti et al.ll20 07: Riekc et al. 2009; Li et alJl2010l ). a good 



morphological correspondence on scales of hundreds of par- 
sees is expected. Fig. [4] indeed shows th at the Ha hot spots 
in the ILR region (Table 4, and also lAlloin et all Il98ll : 
iKristen et ai1ll997r i are clearly associated with IR emitting 
regions (compare with the IRAC 8 /im image as it has the 
best angular resolution). This seems to be the case out to 
A = 100 /xm. At longer wavelengths the decreased angular 
resolution of the Herschel images does not allow us to com- 
pare the morphologies. The opposite is not necessary true 
because the Ha emission is strongly affected by extinction 
caused by the prominent dust lane crossing the ILR region. 
Mid-IR clusters M4 and M5 appear to be detected in the 
Ha image, but they are very faint. This is because they are 
in one of the regions with the highest extinction, as we shall 
see in more detail in Section \4. 2 1 

The Spitzer/IRS map of the [Ne n] 12.81 (j,m emission 
covering the central 2.2 kpc of NGC 1365 is shown in Fig(S] 
This emission line is a good tracer of the current SFR be- 
cause its l uminosity is propo r tional to the numb er of ionizing 
photons l|Roche et al.lll99ll ; lHo~fe Ketoll2007h . The [Nell] 
emission appears as a partial ring and it is enhanced in 
the mid-IR clusters. Its morphology is similar to the PACS 
70 nm morphology and the CO "twin peaks" discussed by 
ISakamoto et~al] \200ih . The AGN of NGC 1365, on the 
other hand is not a bright [Ne u] 12. 81 /xm source, as found 
for ot her active galaxies (see e.g., iPereira- Santacll a et al.l 
12010 j ). Analogously to other SFR indicators, namely the 24 
and 70 /jm emissions, the region to the northeast of the AGN 
(clusters M4, M5, and M6) is brighter in [Nell] than the re- 
gion of clusters M2 and M3. The regions containing these 
bright mid-IR clusters account for approximately 20 — 25% of 
the total [Nell] 12.81 fim emission in the central 2.2 kpc (see 
Table [5}. Note, however that the emission from the clusters 
is unresolved at the Spitzer/IRS SL angular resolution and 
we did not apply a correction for unresolved emission. There- 



2 Note that the Ho+[Nll] image shown in Fig. [4] has not been 
continuum subtracted, and thus we restrict our disc ussion to the 
Ha bright hot spots iden tified by other works (e.g.. lAlloin et al] 
ll98ll;lKristen e t al.lll997t) . 



fore this contribution from the clusters should be taken as 
a lower limit. 

Figure [5] also shows the maps of the 6.2 and the 11.3 /im 
PAH features. Both features probe mostly the emission from 
B stars rather than that from massive on-going star for- 
mation (O stars, see lPeeters et alj|2004 ). The 6.2 /an PAH 
morphology is very similar to [Nell], whereas the 11.3 /im 
PAH map appears to have a deficit of emission in the re- 
gion to the northeast of the AGN. The map of the ratio be- 
tween these two features is in Fig. [7] Although several works 
found variations in the PAH ratios of galaxies not explained 
by changes in the extinction ( see e.g., iGalliano et al.ll2008al ; 
iPereira-Santaella et ahlfeoiObft . in the case of the central re- 
gion of NGC 1365 it might be entirely due to extinction. 
The 6.2 to 11.3 (jm PAH ratio is strongly increased in the 
region enclosing clusters M4, M5, M6, which is also the re- 
gion suffering the highest extinction (see Section 14.21 and 
Fig. 0. Since the relative absorpt ion at 11.3 /jm is h igher 
than at 6.2 /im (see e.g., table 5 in lBrandl et al.| [2006). cor- 
recting the 11.3 /im PAH map for extinction would produce 
a morphology similar to that of the 6.2 /im PAH map. The 
measured EW and PAH ratios of the regions not affected 
by strong extinction are si milar to those measured in other 
starburst galaxies (see e.g., iBrandl et al . 2006; Spoon et al.l 
|2007| ; ISales et alj|2010ft . We conclude that the star forma- 
tion activity inside the ILR region of NGC 1365 has been 
taking place for at least a few tens of million years, as the 
PAH features trace the emission from B stars (jPeeters et al.l 
2004), and is taken place currently because of the bright 
[Nell] 12.81 fim emission requires the presence of young ion- 
izing stars. 



4.2 Extinction: Silicate Feature 

The circumnuclear region of NGC 1365 i s crossed by a 
prominent dust lane entering the ring (see iLindbladl 1 19991 . 
and references therein and also Figs.[T]and[4]|. This dust lane 
is still apparent in the shortest wavelength IRAC bands and 
can be seen passing between clusters M5 and M6 at 3.6 ^m 
(see Fig. [2}. The brightest mid-IR clusters M4, M5, and M6 
are located at the edge of the dust features where the bar 
dust lane enters the ILR region to the northeast of the AGN 
l|Elmegreen et al.|[2009l . and also Figs. |5] and gj. Mid-IR clus- 
ters M2 and M3 to the southwest of the nucleus are in a 
region less affected by extinction and are likely to be associ- 
ated with o ptical super star clu s ters and optically-det ected 
Hn regions (jKristen et al.lll997l ; ISakamoto et al.ll2007h . 

The broad feature at ~ 10 /im is believed to be pro- 
duced by amorphous silicate grains and in normal star- 
form i ng galaxies is observ ed in mild absorption l|Roche et al.l 
Il99ll ; ISmith et alj|200 7bh In the simplest dust geometry of 
a purely absorbing foreground screen the observed apparent 
depth of this feature is proportional to the optical extinc- 
tion. The general variation of the obscuration inside the ILR 
region of NGC 1365 revealed by the optical imaging can be 
traced with the spectral map of the silicate feature. As can 
be seen from Fig. [7J the apparent strength of the silicate 
feature is high (in absorption) in the region encompassing 
clusters M4, M5, and M6 to the northeast of the AGN, al- 
most zero at the location of the AGN, and intermediate (also 
in absorption) around clusters M2 and M3. 

In Table [5] we list the strength of the silicate feature 
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f (100 M m)/f (70 M m) 




Figure 9. The color map is the Herschel/PACS 100 fim to Her- 
schel/PACS 70 fim ratio. The displayed range corresponds to dust 
temperatures of approximately between 22 K (high value of the 
ratio) and 32 K (low value of the ratio). The contours are those 
of the PACS 100 (im image in a square root intensity scale as in 
Fig. [2] The FoV and symbols are also as in Fig. [2] 



measured in the SL ID spectra of selected regi ons. Assum- 
ing a f oreground screen of dust and adopting the lRieke et al.l 
l|l985h extinction law (Ay/Ssi = 16.6), the observed appar- 
ent depths of the star-forming regions in the circumnuclear 
region of NGC 1365 imply values of the visual extinction 
of between Av ~ 4 ± 1 mag and Av ~ 14 ± 1 mag for 
the regions of clusters M2+M3,and M4+M5+M6, respec- 
tively. The high value of the extinction derived for the lat- 
ter region is entirely c onsistent with the values derived by 
iGalliano et ail (|2008ah for the individual clusters using hy- 
drogen recombination lines. Our value of the extinction in 
the region containing M2 and M3 is higher th an the optical 
estima te for the L2 and L3 Hn regions from iKristen et aU 
(I1997D . 

The apparent strength of the silicate feature of the 
AGN measured from the Spitzer/IRS SL spectrum is Ssi = 
—0.10 ±0.05 (Table [SJ) indicating that the feature is present 
slightly in absorption. This value, however, is contaminated 
by extended emission as is apparent from the 8 image 
(see Fig. 3} and the presence of PAH features in the IRS 
nuclear spectrum (Fig. [5]). Indeed, the Gemini/T-ReCS high 
angular resolution spectrum of the AGN is mostly a feature- 
less continuum with no evidence of PAH feature emission 
(Fig. |8]) . Finally, the relatively moderate strength of the sili- 
cate feature of the integrated ~ 30" ~ 2.7 kpc central region 
(S S i = -0.45 ±0.05, Table O of NGC 136^1 is typical of the 
observe d nuclear silicate strengths measured in other local 
LIRG s dPereira-Santaella et al ■ll2010bl ; lAlonso -Herrero et al.l 
12012ft and indicates an average visual extinction in this re- 
gion of Ay ~ 7 ± 1 mag. 



3 The 8 — 1000 ^m IR luminosity of this g alaxy is log(Lm ,/Lpi ) 



4.3 Dust Color Temperatures 

The Herschel images can be used to trace the spatial varia- 
tions of the temperature of the dust inside the ILR region of 
NGC 1365. If we assume that the far-IR emission of a galaxy 
can be approximated with a modified blackbody, then in 
the case of optically thin emission the flux density can be 
expressed as 

f v oc if B(u,T dust ) (2) 

where B(u, Td ust ) is the blackbody function for a dust tem- 
perature of Td us t and /3 is the dust emissivity. In the simplest 
approximation we can calculate the color temperature T c of 
the dust using the ratio of the surface brightness at two 
wavelengths with the following equation expressed in terms 
of wavelengths 



3+/3 



3 hc/A 2 KT c 
>hc/\ 1 KT c 



(3) 



11.03 using the IRAS flux densities from Sanders et 



al.l jgOOgj) . 



We chose to construct a map of the PACS 100 fim to the 
PACS 70 fim ratio as these two wavelengths provide the best 
angular resolutions with Herschel. We rebinned the PACS 
70 fim image to the pixel size of the PACS 100 fim image. 
We matched the PSF of the two images by smoothing the 
70 fim with a Gaussian funct ion, although this may intro- 
duce some artifacts (see e.g., iBendo et al. I l201Cl 12012ft . We 
then calculated the dust color temperature by solving Equa- 
tion |3] We fixed the value of the dust emissivity /3 — 2, as it 
is found to fit well the integrated SEDs of local IR galaxies 
l|Dunne fc Eales!l200ll ). 

The map of the PACS 100 fim to 70 fim ratio of 
the inner ~ 60" is shown in Fig [5] The observed 
range of /„(100 fim) / f ^(70 fim) in the central region trans- 
lates into values of the dust color temperature of be- 
tween T c (100/im/70/im) ~ 32 K and T c (100/im/70/im) ~ 
22 K. Within the ILR region of NGC 1365 the highest 
T c (100/im/70/im) color temperatures correspond to regions 
actively forming stars, that is, the regions containing the 
mid-IR star clusters and the L4 Ha hot spot. These bright 
star forming regions have color temperatures similar to, al- 
though slightly higher than, those of H n regions in the spiral 
arms. The regions with the coldest T c (100fim/70fim) in the 
ILR region might be associated with some of the foreground 
dust features seen in the optical images (see Fig.[l|. Finally, 
the region of the AGN does not appear different in terms 
of the color temperature when compared to the bright star 
forming regions. This is probably because the AGN is faint 
in the far-IR (see next section) coupled with the relatively 
coarse angular resolution of this map (i.e., that of the PACS 
100 /xm). 



4.4 Dust Properties 

In this section we study the properties of the dust in the 
ILR region of NGC 1365, particularly the dust heated by 
processes other than the AGN. To do so, we first subtracted 
the AGN emission (see Table [3]) from the observed fluxes in 
the Spitzer/MIPS 24 fim band and all the Herschel bands. 
We then fitted the observed non-AGN mid- and far-IR SED 
using a combination of t wo modified blackbody following 
iDunne fc Ealesl (120011 ) and lClements et all (|2010ft : 

/„ = B(v,T„) + N c v p B{v,T c ) (4) 
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Figure 10. Best fit (solid line) to the AGN-subtracted SED of 
the IR emission in the ILR region of NGC 1365 (r = 30", filled 
dots). The best fit was achieved using two modified blackbodies 
with temperatures T c = 24 K and T w = 54 K (dashed lines) . The 
dotted lines are fits done by varying T c by ±1 K, and keeping the 
temperature of the warm dust fixed. 



where iV w and N c are the relative masses of the warm and 
cold dust components and T c and T w are their temperatures. 
As done in Section 14.31 we fixed the dust emissivity /3 = 2 
and then used a standard \ 2 minimization technique to fit 
the two dust temperatures and the relative masses of the 
two dust components. For the fit we used the fluxes mea- 
sured through an r = 30" aperture, and normalized them to 
that at 100 /mi. For the Spitzer/MlPS 24 /im flux density, 
we added in quadrature the photometric error and the un- 
certainty associated with subtracting the AGN component 
at this wavelength to a total error budget of ~ 20%. For the 
Herschel data points because the AGN emission in the far- 
IR is very small compared to the total emission inside the 
ILR region, the dominant source of error is that associated 
with the photometric calibration of the data: ~ 10%. 

Fig. UPl shows the best fit to the AGN-subtracted SED 
of the ILR region of NGC 1365, which was obtained with 
T c = 24 K and T w = 54 K, and relative masses of cold to 
warm dust of N c /Nw ~ 120. It is worth noting that the 
fit using two modified blackbodies is formally better than a 
fit to the far-IR data (i.e., excluding the 24 fj,m data point) 
using a single modified blackbody. This can be explained be- 
cause there is a non-negligible contribution from the warm 
dust component to the Herschel PACS 70 fjxa flux density, 
as found for t he integrated emis s ion of other nearby galax- 
ies (see e.g., iBendo et ail |2010| ; ISmith et all l2O10h . How- 
ever, the temperature of the warm component is not tightly 
constrained since varying the warm dust temperature by 
as much as ±4K produces statistically similar good fits 
(x 2 ^ 2 Xmin)- This is because the peak of this component at 
around 50 /im is not well sampled with the present data and 
the mass contribution of the warm component is small. The 
temperature of the cold component, on the other hand is well 
constrained. To get an estimate of the uncertainties on the 
fitted cold dust temperature, we changed it while fixing the 
temperature of the warm component. As can be seen from 
Fig. 1101 cold dust temperatures in the range T c = 24 ± 1 K 
provide acceptable fits to the SED. 



Previous works showed that the integrated mid-IR and 
far-IR SEDs of Seyfert galaxies could be fitted with a 
combination of various dust temperatures probing differ- 
ent heating mechanisms. These include warm dust associ- 
ated with the AGN torus, cold dust similar to that observed 
in starburst galaxies, and very cold dust at temperatures 
characteristic of dust heated by the inte r stellar medium 
llPerez Garcia fc Rodriguez Espinosal l200ll ; ISpinoglio et~ail 
2002). With the superior angular resolution of the Herschel 
observatory we can now obtain spatially resolved observa- 
tions of nearby Seyfert galaxies and star forming galax- 
ies. The dust temperature of the ring of star formation in 
NGC 1365 is similar to that of other nuclear and circum- 
nuclear starbursts with or without an AGN with similar 
IR luminos ities and physical sizes, for i nstance, the ring of 
NGC 3081 dRamos Almeida et al.l2011al ) and the nuclear re- 
gion of M83 (|Fovle et alj|2012h . On the other hand, the dust 
of Mrk 938, which is part of our survey of Seyfert galaxies, 
has a considerably higher dust temperature of T = 36±4K, 
probably due to the smaller size of the IR emitting region of 
this galaxy and higher IR luminosity (Esq uei et al.ll2012i ). 

We calculated the dust mass (Md ust ) within the ILR 
region of NGC 1365 using the foll owing equation fr om 
IClements etaH (|201Ch (adapted from iHildebrand II1983T ): 



M dust = 



N c 



+ ■ 



fvD 2 

Kdust(^) " \B(u,T c ) ' B(u,T v ) 



(5) 



where /„ is the observed flux density, D is the luminosity 
distance, B(v,T) is the blackbody emission for the best fit- 
ting dust temperatures, and Kd(^) i s the d ust absorption 
coefficient. As done by lEsquei et al.l ([2012) , we evaluated 
this expression at 250 ^m using an absorption coefficient of 
Kdust(250 um) = 4.99cm 2 g~\ as i nterpolated from the dust 
model of (Li fc Draine 200 ll l|200ll ). We derived a dust mass 
in the ILR region of M dust (ILR) = 6.9 x 10 r M , which ac- 
counts for approximately 25% of the total dust mass of this 
gala xy Md U3 t (total) = 3 x 10 s Mq. The latter estimate is 
from lWiebe et ail l|2009l ) but recalculated for the absorption 
coefficient used in this work. 



4.5 Star Formation Rate 

The ages an d masses of the mid-IR clusters (~ 6 — 8Myr and 
~ 1O 7 M0, iGalliano et al.|[2008al ) and the large amount of 
molecular gas available in the central regions of this galaxy 
l|Sakamoto et al.l 120071 ) indicate that there is intense on- 
going star formation activity there. We can use the IR ob- 
servations of NGC 1365 to estimate the obscured SFR in the 
ILR region, and compare it with th e unobscured SFR. As 
discussed by iKennicutt et al.l (|2009h , a combination of the 
observed (not corrected for extinction) Ha luminosity and 
a mid-IR monochromatic luminosity (preferably 24 fim) will 
provide the best estimate of the total SFR in a moderately 
obscured environment (see Section \A.2\ , such as the ILR re- 
gion of NGC 1365. We note that this empirically calibrated 
recipe includes contributions fro m dust heating from al l stars 
and not only the youngest (see IKennicutt et al.l 2009|. for a 
full di scussion). Using the Ha flux of Forster Schreiber et al.l 
l|2004h for a 40"-diameter aperture and our 24 /im flux (af- 
ter subtracting the AGN component) we estimated a to- 
tal SFR within the ILR region (inner ~ 5.4 kpc, r = 30") 
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of NGC 1365 of SFR= 7.3M yr _1 for a Salpeter initial 
mass function (IMF). Approximately 85% of this SFR is 
contributed by the 24 /xm emission, and thus it originates 
from dust-obscured star forming regions. 

Most of the IR emission in the ILR region comes from 
the star formati on ring contain i ng th e bright mid-IR clus- 
ters identified bv lGalliano et all (|2005l ). Indeed, we estimate 
that within the inner 2.7kpc the total SFR is 5.6M0yr _1 . 
This has been calculated using the AGN-subtracted 24 fim 
flux density ( r = 15") and the 24 " -diam eter aperture 
Ha flux from iForster Schreiber et al.l |2004l ). This corre- 
sponds to a SFR surface density in the ring of Esfr = 
2.2 Mq yr _1 kpc -2 , for a ring radius of 900 pc. This value of 
the SFR d ensity is similar t o those of other circumnuclear 
starbursts jKennicutt Ill998l ) and is expected to be larger 
by factors of 100 — 1000 compared to the disk of the galaxy 
l|Elmegreen II 1994! ) . 



5 SUMMARY AND CONCLUSIONS 

In this paper we have studied the IR emission associ- 
ated with the star formation activity in the ILR region 
of NGC 1365, as well as the IR emission of the AGN. To 
this end we have analyzed new far-IR (70 — 500 /im) Her- 
schel/PACS and SPIRE imaging, and high angular resolu- 
tion (~ 0.4") Gemini/T-ReCS mid-IR imaging and spec- 
troscopy of this galaxy. We have also made use of archival 
Spitzer/IRAC and MIPS imaging and IRS spectral mapping 
data. Our main findings for the inner D ~ 5 kpc region of 
NGC 1365 are: 

• The new Herschel PACS imaging data at 70, 100, and 
160 /im reveal that the ring of star formation in the ILR 
region is bright in the far-IR. The AGN is the brightest 
mid-IR source in the inner 2 kpc up to A ~ 24 /im, but it 
becomes increasingly fainter in the far-IR when compared 
with the mid-IR clusters or groups of them in the ring. 

• The 24 and 70 /im emissions as well as the 
[Nell] 12.81 fim line and PAH features trace the star-forming 
ring in the ILR region and have morphologies similar to the 
CO "twin-peaks". This all indicates that there is intense 
on-going star formation taking place in the inner few kpc of 
NGC 1365. 

• The unresolved near and mid-IR nuclear emission 
and mid-IR spectrum (i.e., AGN-dominated emission) of 
NGC 1365 are well reproduced with a relatively compact 
torus (outer radius of R = 5^ 5 P c ) with an opening an- 
gle of atoms = 36lg 4 deg, and an AGN bolometric luminosity 
L bo i(AGN) = 2.6±0.5xl0 43 ergs -1 using the clumpy torus 
models. These parameters are in good agreement with inde- 
pendent estimates in the literature. 

• Using the fitted torus model we quantified the AGN 
emission in the far-IR. The AGN only contributes at most 
1% of the 70 fim emission within the inner 5.4 kpc (r — 30"), 
and less than 1% at longer wavelengths. At 24 /im the AGN 
accounts for ~ 15% of the emission in the same region. We 
estimated that the AGN bolometric contribution to the 3 — 
1000 /im luminosity in the inner 5.4 kpc is approximately 
5%. 

• The non-AGN 24 to 500 /im SED of the ILR region (in- 
ner 5.4 kpc) of NGC 1365 is well fitted with a combination of 
two modified blackbodies with warm and cold temperatures 



of 54 K and 24 K, respectively. However, the cold dust com- 
ponent accounts for most of total dust mass inferred in this 
region (Md us t(ILR) = 7 x 10 7 Mq) and has a temperature 
similar to that of other nuclear and circumnuclear starbursts 
of similar sizes and IR luminosities. 

• From the comparison between the SFR from Ha (un- 
obscured) and the SFR from 24 pirn (obscured) we infer that 
up to ~ 85% of the on-going SFR inside the ILR region of 
NGC 1365 is taking place in dust-obscured regions in the 
ring of star formation. 
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